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O ! Abstract 



The design of Medium Access Control (MAC) protocols for wireless sensor networks (WSNs) has 

been conventionally tackled by assuming battery-powered devices and by adopting the network lifetime 

as the main performance criterion. While WSNs operated by energy-harvesting (EH) devices are not 

. limited by network lifetime, they pose new design challenges due to the uncertain amount of harvestable 

energy. Novel design criteria are thus required to capture the trade-offs between the potentially infinite 

' network lifetime and the uncertain energy availability. 

ON ■ 



This paper addresses the analysis and design of WSNs with EH devices by focusing on conventional 
MAC protocols, namely TDMA, Framed-ALOHA (FA) and Dynamic-FA (DFA), and by accounting 
for the performance trade-offs and design issues arising due to EH. A novel metric, referred to as 
delivery probability, is introduced to measure the capability of a MAC protocol to deliver the measure 
, of any sensor in the network to the intended destination (or fusion center, FC). The interplay between 

^ , delivery efficiency and time efficiency (i.e., the data collection rate at the FC), is investigated analytically 

using Markov models. Numerical results validate the analysis and emphasize the critical importance of 
accounting for both delivery probability and time efficiency in the design of EH- WSNs. 
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I. Introduction 

Recent advances in low-power electronics and energy-harvesting (EH) technologies enable 
the design of self-sustained devices that collect part, or all, of the needed energy from the 
surrounding environment. Several systems can take advantage of EH technologies, ranging from 
portable devices to wireless sensor networks (WSNs) |[T1. However, EH devices open new design 
issues that are different from conventional battery-powered (BP) systems [2], where the main 
concern is the network lifetime In fact, EH potentially allows for perpetual operation of the 
network, but it might not guarantee short-term activities due to temporary energy shortages [[2l. 
This calls for the development of energy management techniques tailored to the EH dynamics. 
While such techniques have been mostly studied at a single-device level flU, in wireless scenarios 
where multiple EH devices interact with each other, the design of EH-aware solutions needs to 
account for a system-level approach [I5]|[l6]|. This is the motivation of this work. 

In this paper, we focus on system-level design considerations for WSNs operated by EH- 
capable devices. In particular, we address the analysis and design of medium access control 
(MAC) protocols for single-hop WSNs (see Fig. [1]) where a fusion center (EC) collects data 
from sensors in its surrounding. Specifically, we investigate how performance and design of MAC 
protocols routinely used in WSNs, such as TDMA [7J, Eramed- ALOHA (EA) and Dynamic-EA 
(DEA) (HI, are influenced by the discontinuous energy availability in EH-powered devices. 

A. State of the Art 

In recent years, WSNs with EH-capable nodes have been attracting a lot of attention, also 
at commercial level. To provide some examples, the Enocean Alliance proposes to use a MAC 
protocol for EH devices based on pure ALOHA strategies [[TT|. while an enhanced self-powered 
REID tag created by Intel, referred to as WISP [[T2|. has been conceived to work with the EPC 
Gen 2 standard [T31 that adopts a EA-like MAC protocol. 

However, while performance analysis of MAC protocols in BP-WSNs have been investigated 
in depth (see e.g., [|7l[[8l[[T4|). analyses of MAC protocols with EH are hardly available. A notable 
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Figure 1. WSN with a single Fusion Center (FC) gatliering data from M sensors, whicli are equipped witli an energy storage 
device (BSD) and an energy-liarvesting unit (EHU). 

exception is ["61, where data queue stability has been studied for TDMA and carrier sense multiple 
access (CSMA) protocols in EH networks. We remark that routing for EH networks has instead 
received more attention, see e.g., ||2l lfT5l . 

B. Contributions 

In this paper we consider the design and analysis of TDMA, FA and DFA MAC protocols 
in the light of the novel challenges introduced by EH. In Sec. Un] we propose to measure the 
system performance in terms of the trade-off between the delivery probability, which accounts 
for the number of sensors' measurements successfully reported to the FC, and the time efficiency, 
which measures the rate of data collection at the FC (formal definitions are in Sec. Un]). We then 
introduce an analytical framework in Sec. |IV]and Sec. |V]to assess the performance of EH-WSNs 
in terms of the mentioned trade-off for TDMA, FA and DFA MAC protocols. In Sec. |Vl] we 
tackle the critical issue in ALOHA-based MAC protocols of estimating the number of EH sensors 
involved in transmission, referred to as backlog, by proposing a practical reduced-complexity 
algorithm. Finally, we present extensive numerical simulations in Sec. IVIII to get insights into 
the MAC protocol design trade-offs, and to validate the analytical derivations. 
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II. System Model 

In this paper, we consider a single-hop WSN with a FC surrounded by M wireless sensors 
labeled as S*!, S'2, Sm (see Fig. [I])- Each sensor (or user) is equipped with an EH unit (EHU) 
and an energy storage device (ESD), where the latter is used to store the energy harvested by the 
EHU. The FC retrieves measurements from sensors via periodic inventory rounds (IRs), once 
every Tint seconds [s\. Each IR is started by the FC by transmitting an initial query command (Q), 
which provides both synchronization and instructions to sensors on how to access the channel. 
Time is slotted, with each slot lasting [s] . The effective duration of the nih. IR, during which 
communication between the FC and the sensors takes place, is denoted by TiR{n). We assume 
that Tjfj{n) <^ Tint for all IR n, and also that the query duration is negligible, so that the ratio 
Tjfj{n)/Ts indicates the total number of slots allocated by the FC during the nth IR. 

In every IR, each sensor has a new measure to transmit with probability a, independently 
of other sensors and previous IRs. If a new measure is available, the sensor will mandatory 
attempt to report it successfully to the FC as long as enough energy is stored in its ESD (see 
Sec. III-BI for details). Each measure is the payload of a packet, whose transmission fits within 
the slot duration Ts. Sensors' transmissions within each IR are organized into frames, each of 
which is composed of a number of slots that is selected by the FC. Depending on the adopted 
MAC protocol, any user that needs to (and can) transmit in a frame either chooses or is assigned 
a single slot within the frame for transmission as it will be detailed below. Moreover, after a 
user has successfully transmitted its packet to the FC, it first receives an acknowledge (ACK) 
of negligible duration by the FC and then it becomes inactive for the remaining of the IR. We 
emphasize that the FC knows neither the number of sensors with a new measure to transmit, 
nor the state of sensors' ESDs. 



A. Interference Model 

We consider interference-limited communication scenarios where the downlink packets trans- 
mitted by the FC are always correctly received (error-free) by the sensors, while uplink packets 
transmitted by the sensors to the FC are subject to communication errors due to possible 
interference arising from collisions with other transmitting sensors. The uplink channel power 
gain for the mth sensor during the nth IR is hm{n). Channel gain hm{n) is assumed to be 
constant over the entire IR but subject to random independent and identically distributed (i.i.d.) 
fading across IRs and sensors, with pdf fh{-) and normalized such that E [hm{n)] = 1, for all 
n, m. In the presence of simultaneous transmissions within the same slot during the kth frame 
of the nth IR, a sensor, say S^, is correctly received by the FC if and only if its instantaneous 
signal-to-interference ratio (SIR) 7^,^ (n) is larger than a given threshold 'jth, i-e-^ if 

lm,k [n) = ^ T-r^ > 1th, (1) 

where Xm,k{n) denotes the set of sensors that transmit in the same slot selected by Sm in frame 
k and IR n. We assume 7^/1 > QdB so that, in case a slot is selected by more than one sensor, 
at most one of the colliding sensor can be successfully decoded in the slot. 

According to the interference model ([T]), any slot can be: empty when it is not selected by 
any sensor; collided when it is chosen by more than one sensors but none of them transmits 
successfully; successful when one sensor transmits successfully possibly in the presence of other 
(interfering) users. Successful transmission in the presence of interfering users within the same 
slot is often referred to as capture effect lfT4l . 

Remark 1: Errors in the decoding of downlink query packets can be accounted for through 
the parameter a as well. In fact, let aq be the probability that a user correctly decodes the 
downlink packet sent by the FC at the beginning of an IR. Moreover, assume that downlink 
decoding errors are i.i.d. across sensors and IRs, and let a; at be the probability that a user has a 
new measure to transmit in any IR. Then, the probability that any user Sm has a new measure 
and correctly decodes the FC's query is given by the product a = aqa^- 
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B. BSD and Energy Consumption Models 

We consider a discrete ESD with + 1 energy levels in the set £ = {0, S, 25, N5}, where 
5 is referred to as energy unit. Let Em{n) E S he the energy stored in the ESD of the mth 
user at the beginning of the nth IR. Energy Em{n) is a random variable that is the result of the 
EH process and the energy consumption of the sensor across IRs; its probability mass function 
(pmf) is pE(n) (■) and the corresponding complementary cumulative distribution function (ccdf) 
is GE(n) (x) = Fi[Em{n) > x]. Note that, the initial energy distribution pe(i) (•) is given, while 
the evolution of the pmf pE{n) (■) for n > 1 depends on both the MAC protocol and EH process. 

We assume that each time a sensor transmits a packet it consumes an energy e, which accounts 
for the energy consumed in the: a) reception of the EC's query that starts the frame (see Eig. 
[21); b) transmission; c) reception of EC's ACK or not ACK (NACK) packet. At the beginning of 
each IR, a sensor with a new measure to transmit can participate to the current IR only if the 
energy stored in its ESD is at least e. Let es = e/S he the number of energy units S required 
for transmission, where ss is assumed to be an integer value without loss of generality. Let 
= N5/e = N/es he the (normalized) capacity of the ESD, which is assumed to be an integer 
indicating the maximum number of (re)transmissions allowed by a fully charged ESD. 
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Figure 2. Examples for TDMA and DFA MAC protocols for Af — 3. FA is not depicted since it is a special case of DFA 
with only one frame. The backlog for each frame is indicated above each query. Some sensors might not be in the backlog due 
to energy shortage and/or absence of a new measure to report. 



C. Energy Harvesting Model 

During the time T^t between the nth and (n + l)th IRs the mth sensor 5*^ harvests an energy 
EH,m{n), which is modeled as a discrete random variable, i.i.d. over IRs and sensors, with pmf 
Qi = FilEH.min) = i6], with i E {0, 1, 2, ...}, and for all m and n. For technical reasons that we 
discuss in Sec. IV-B[ we assume that the probability go and qi of harvesting zero and one energy 
unit respectively, are both strictly positive, namely go > and gi > 0. 

We assume that the EH dynamics is much slower than the IR duration TjR(n), so that the 
amount of energy harvested within Tjji(n) can be considered as negligible with respect to e (recall 
also that Tjji{n) <^ Tint)- Hence, the only energy that a sensor can actually use throughout an 
IR is the energy initially available at the beginning of the IR itself (i.e., Em{n)). 



III. Performance Metrics and Medium Access Control Protocols 

We first introduce in Sec. IIII-AI the considered performance metrics, namely delivery proba- 
bility and time efficiency, and then in Sec. IIII-BI we review the considered MAC protocols. 
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A. MAC Performance Metrics 

1 ) Delivery Probability: The delivery probability (n) measures the capability of the MAC 
protocol to successfully deliver the measure of any sensor, say Sm, to the FC during the nth IR 

Pd{n) = Pr [Sm transmits successfully in IR n| has a new measure in IR n] . (2) 

The statistical equivalence of all sensors makes the probability Q independent of the specific 
sensor. Notice that a sensor fails to report its measure during an IR if either it has an energy 
shortage before (re)transmitting the packet correctly, or the MAC protocol does not provide 
the sensor with sufficient retransmission opportunities. Given the potentially perpetual operation 
enabled by EH, it is relevant to evaluate the delivery probability when the system is in steady- 
state. The asymptotic delivery probability is thus obtained by taking the limit of pd (n) for large 
IR index n, provided that it exists, as 

pf = limpd(n). (3) 

2) Time Efficiency: The time efficiency pt{n) measures the probability that any slot allocated 
by the MAC within the nth IR is successfully used (i.e., it is neither empty nor collided, see 
SecHEAl) 

Pt (n) = Pr [The FC correctly retrieves a packet in any slot of the nth IR] . (4) 

By taking the limit of (51) for n — )■ oo, we obtain the asymptotic time efficiency 

pf^ = lim Pt (n) . (5) 

Remark 2: Informally speaking, the time efficiency Pt{n) measures the ratio between the total 
number of packets successfully received by the FC and the total number of slots allocated by the 
MAC protocol (i.e., TjR{n)/Ts, see Sec. HB- As it will be shown in Sec. IIII-B[ the IR duration 
Tmin) is in general a random variable, and consequently, time efficiency pt{n) differs from 
more conventional definitions of throughput (see e.g., [[8]|) which measure the number of packets 
delivered over the interval between two successive IRs Tint, instead of Tjnin). The rationale 
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for this definition of time efficiency is that it actually captures more effectively the rate of data 
collection at the FC. Whereas, the delivery probability accounts for the fraction of users, with 
a new measure to transmit at the beginning of the current IR, which are able to successfully 
report their payload to the FC within the IR, where delivery failures are due to collisions and 
energy shortages. 

In contention based MACs (e.g., ALOHA), there is a trade-off between delivery probability 
and time efficiency. In fact, increasing the former generally requires the FC to allocate a larger 
number of slots in an IR to reduce packet collisions, which in turn decreases the time efficiency. 

B. MAC Protocols 

In this section, we review the standard MAC protocols that we focus on. 

1) TDMA: With the TDMA protocol, each user is pre-assigned an exclusive slot that it 
can use in every IR, irrespective of whether it has a measure to deliver or enough energy to 
transmit. Recall that such information is indeed not available at the FC. Every nth IR is thus 
composed by one frame with M slots and has fixed duration T™ = MTs, as shown in Fig. [21 
Since TDMA is free of communication errors in the considered interference-limited scenario, its 
delivery probability Pdin) is only limited by energy availability and it is thus an upper bound 
for ALOHA-based MACs. However, TDMA might not be time efficient due to the many empty 
slots when the probability of having a new measure a and/or the EH rate are small. 

2) Framed-ALOHA (FA) and Dynamic-FA (DFA): Hereafter we describe the DFA protocol 
only, since FA follows as a special case of DFA with no retransmissions capabilities as discussed 
below. The nth IR, of duration Tfjf'^{n), is organized into a set of frames as shown in Fig. 
[2l The backlog Bk{n) for the A;th frame is the set composed of all sensors that simultaneously 
satisfy the following three conditions: z) have a new measure to transmit in the nth IR; ii) have 
transmitted unsuccessfully (because of collisions) in the previous k — \ frames (this condition 
does not apply for frame k = 1); Hi) have enough energy left in the ESD to transmit in the 
/cth frame. All the users in the set Bk{n), whose cardinality |-BA;(n,)| = Bk{n) is referred to as 
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backlog size, thus attempt transmission during frame k. To make this possible, the FC allocates 
a frame of Lk{n) slots, where Lk{n) is selected based on the estimate Bk{n) of the backlog size 
Bk{n) (estimation of Bk{n) is discussed in Sec. IVTl) as 



Lk(n) 



pBkin) 



(6) 



where [■] is the upper nearest integer operator, and p is a design parameter. Note that, if the 
backlog size is B, the probability /3 (j, B, L) that j < B sensors transmit in the same slot in a 
frame of length L is binomial [fT6l 



Finally, FA is a special case of DFA where only one single frame of size Li (n) is announced 
as retransmissions are not allowed within the same IR. 

IV. Analysis of the MAC Performance Metrics 

In this section we derive the performance metrics defined in Sec. IIII-AI for TDMA, FA and 
DFA. The analysis is based on two simplifying assumptions: 

A.l Known backlog: the FC knows the backlog size Bk{n) = \Bk{n)\ before each fcth frame; 

A.2 Large backlog: the backlog size Bk{n), in any IR n and any frame k of size Lk{n) = 
\pBk{n)~\, is large enough to let the probability dV]) be approximated by the Poisson distri- 
bution 1161: 

P{j,Bkin),L,in))^^. (8) 

Assumption A.l simplifies the analysis as in reality the backlog can only be estimated by the 
FC (see Sec.|Vl]and Sec. I VIII for the impact of backlog estimation). Assumption A.2 is standard 
and analytically convenient, as it makes the probability (3 (j, Bk{n), Lk{n)) dependent only on 
the ratio p between the frame length Lk{n) and the backlog size Bk{n). The assumptions above 
are validated numerically in Sec. IVIII 
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A. Delivery Probabilities 

Here we derive the delivery probability ^ within any nth IR under the assumptions A.l and 
A.2 for the considered MAC protocols. The IR index n is dropped to simplify the notation. 

1) Delivery Probability for TDMA: As the TDMA protocol is free of collisions, each sensor 
Sm that has a new measure to report in the current IR cannot deliver its payload to the FC only 
when it is in energy shortage, namely if Em < s. Provided that user Sm has a new measure to 
transmit, the delivery probability ^ reduces to 

= Pr [Em >e] = {e) , (9) 

which is independent of the sensor index m and dependent only on the ccdf (■) of the 
energy stored in sensor BSD at the beginning of the considered IR. The BSD energy distribution 
for any arbitrary nth IR is derived in Sec. |Vl 

2) Delivery Probability for FA: In the FA protocol, each sensor Sm that has a new measure 
to report in the current IR is able to correctly deliver its payload to the FC only if: a) it transmits 
successfully in the selected slot, possibly in the presence of interfering users provided that its 
SIR is 7m,i > 7t/i; and b) it has enough energy to transmit. From ©, the probability that sensor 
Sm, with Sm G Bi, transmits successfully in the selected slot, given that |Xm,i| = j users select 
the same slot of Sm (thus colliding), is given by 



Pc{]) = Pr 



(10) 



hm > 1th ^ hi 

1=1 

where, without loss of generality, we assumed that Xm,i = {Si, Sj}, and Sm ^ 2^1, as users 
are stochastically equivalent. Under the large backlog assumption A.2, the probability that there 
are j interfering users is Poisson-distributed (see ([8])), and thus the unconditional probability pc 
that Sm captures the selected slot can be approximated as 

p,~e-^^--p,(j). (11) 
j=o ^ 

Note that, in (fTT)) we also extended the number of possible interfering users up to infinity as 
Pc{j) rapidly vanishes for increasing j. Moreover, depending on the channel gain pdf fh{-). 
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probabilities (flOl) can be calculated either analytically (e.g., when //,(■) is exponential, see ifTTll ) 
or numerically. 

Finally, under assumption A.2, the successful transmission event is independent of the BSD 
energy levels (which in principle determine the actual backlog size in (|7])), and thus the delivery 
probability (O for the FA protocol can be calculated as the product between the probability 
(^) = l^m > that sensor Sm has enough energy to transmit and the (approximated) 
capture probability (fTTI) as 

oo _ 

where the ESD energy ccdf G^^ (e) for any arbitrary nth IR is derived in Sec. |Vl 

3) Delivery Probability for DFA: DFA is composed of several instances of FA, one for each 
kth frame of the current IR. As DFA allows retransmissions, we need to calculate the probability 
Pc,k{j) that any sensor active during frame k, say Sm £ ^k, transmits successfully in the selected 
slot given that there are \lm,k\ = j users that transmit in the same slot, with Xm,k ^ -Bfc. The 
computation of Pc,k{j), for A; > 1, is more involved than (flOl) . In fact, packets collisions introduce 
correlation among the channel gains of collided users, as any sensor in the backlog S^, for /c > 1, 
might have collided with some other sensors in the set Bk- We recall that, even though the channel 
gains are i.i.d. at the beginning of the IR, they remain fixed for the entire IR. 

While the exact computation of probabilities Pc,k{j) is generally cumbersome, the large backlog 
assumption A.2 enables some simplifications. Specifically, correlation among channel gains can 
be neglected, since for large backlogs it is unlikely that two users collide more than once within 
the same IR. By assuming independence among the channel gains at any frame, calculation of 
Pc,k{j) requires only to evaluate the channel gain pdf fh^\-) at the A;th frame for any user within 
Bk, which is the same for all users by symmetry. The computation of pdf fl''\-) can be done 
recursively, starting from frame A; = 1, so that at frame k we condition on the event that the 
SIR m was 7m,fc-i < 1th- Under assumption A.2, this can be done numerically. 

Now, let hm\ for m G {1, ...,M} and k G {1, ...,F^}, be random variables with pdf fl^\-) 
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independent over m, where hm = hm- The conditional capture probabilities Pc,k{j) can then be 
approximated as (compare to (flOl) ) 

(13) 

for any m ^ as users are stochastically equivalent. By exploiting the Poisson ap- 

proximation similarly to (fTTI) . the unconditional probability that any user within the backlog 
successfully transmits in the selected slot during the kth frame becomes 

oo ^ 

Pc,k ^ e~p V —rrPckU)- (14) 
j=o ^ •> 

Recalling that a user keeps retransmitting its message until it is successfully delivered to the 
FC, then the successful delivery of a message in a frame is a mutually exclusive event with respect 
to the delivery in previous frames. Therefore, the probability of transmitting successfully in the 
/cth frame, given that enough energy is available, is pc,fc Ilti (1 ~ Vc,i) ■ Finally, by accounting 
for the probability {ke) = Pr [Em > ke] of having enough energy in each fcth frame, the 

DFA delivery probability can be obtained, under assumption A.2, a^l 

k-l 

^ J2 G^/^ ike) p,,, n (1 - P..) > (15) 

fc=l i=\ 

where the BSD energy ccdf G^^^ (ke) for any arbitrary nth IR is derived in Sec. |Vl 
B. Time Efficiencies 

In this section we derive the time efficiency dH) for the three considered protocols. 

1) Time Efficiency for TDMA: Let be the event indicating that user SjYi has a new 
measure to report in the current IR, with Pr[A^m] = en, then the TDMA time efficiency (HI) is 
given by the probability that the mth user has enough energy to transmit and a new measure to 
report: 

= Pr [Em > e, Mm] = Pr [Em > e] Pr [Mm] = aCl"" (e) , (16) 
where we exploited independence between energy availability Em and Mm- 

'Note that in principle the backlogs Bi,B2--- are correlated, and therefore the exact p^^^ should be obtained by averaging 
over the joint distribution of the backlog sizes. However, the assumption A.2 removes the dependence on the backlog size. 
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2) Time Efficiency for FA: Since we assumed 'yth > OdB, then when more than one user 
transmits within the same slot, only one of them can be decoded successfully, that is, successful 
transmissions of different users within the same slot are disjoint events. Therefore, the probability 
that a slot, simultaneously selected by j users, is successfully used by any of them is given by 
jPc{j — 1)5 where pc{j — 1) is (flOl) by recalling that any user have (j — 1) interfering users. 
Furthermore, under assumption A.2, the probability that exactly j users select the same slot is 

e~p/ [fP and by summing up over the number of simultaneously transmitting users j we get 

00 00 

j=i -J- j=o 
Note that, a consequence of assumption A.2 is to make the FA time efficiency (flTl) independent 

of the ESD energy distribution. Moreover we remark that, when p = 1, pc{j) = 1 for j = and 

PcU) = for j > 0, then we have pf^ = e^^, which is the throughput of slotted ALOHA dHl. 

3) Time Efficiency for DFA: The derivation of the DFA time efficiency pf^^ follows from 
the FA time efficiency by accounting for the presence of multiple frames within an IR similarly 
to Sec. IIV-A3[ Since the time efficiency is defined over multiple frames, we first derive the time 

efficiency in the /cth frame, similarly to (fT71) but considering (fT3l) instead of (flOl) . as 

00 ^ 



We then calculate pf^"^ by summing (fTSi) up, for all k E {l,...,Fs}, weighted by the (ran- 
dom) length of the corresponding frame normalized to the total number of slots in the IR 
'^k=i ^k- Note that, under assumption A.2 the random frame length is well-represented by 
its (deterministic) average value Lk c:^ E [Lk] = pE [Bk\ and thus the DFA time efficiency results 



where the average backlog size E[Bk\ in frame k, can be computed, under assumption A.2, as 
E[Bk] = MaG'g^^(fe) nti^ (1 - Pc,^)■ In fact. Ma indicates the average number of users that 
have a new measure to report in the current IR, G{ke) is the probability that ke energy units 
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are stored in the ESD at the beginning of the IR, thus allowing k successful transmissions, and 
Y[i=i (1 ~ Pc,i) is the probability that a sensor collides in all of the first (A; — 1) frames. 

V. ESD ENERGY EVOLUTION 

In Sec. |IV] we have shown that the performance metrics for the nth IR depend on the energy 
distribution in the sensor ESD at the beginning of the IR itself. The goal of this section is to 
derive the ccdf GE{n){-), for any IR n, in order to obtain the asymptotic performance metrics 
© and © from Sec. IIV-AI and Sec. |I\^ respectively. 

In general, the evolution of sensor ESDs across IRs in DFA are correlated with each other, due 
to the possibility of retransmitting after collisions. However, under the large backlog assumption 
A.2, similarly to the discussion in Sec. IIV-A3[ the evolution of sensor ESDs become decoupled 
and can thus be studied separately. Accordingly, we develop a stochastic model, based on a 
discrete Markov chain (DMC) that focuses on a single sensor ESD as shown in Fig. [3l In 
addition, we concentrate on the DFA protocol as ESD evolutions for TDMA and FA follow as 
special cases. Note that, in TDMA (or FA), the evolution of sensor ESDs are actually independent 
with each other as retransmissions are not required (or allowed). 

A. States of a Sensor 

The state of a sensor is uniquely characterized by: /) sensor activity or idleness (see below); 
//) the amount of energy stored in its ESD; ///) the current frame index if the sensor is active. 
A sensor is active if it has a new measure still to be delivered to the EC in the current IR and 
enough energy in its ESD, while it is idle otherwise. States in which the sensor is active, referred 
to as active states, are denoted by and they are characterized by: a) the current frame index 
k G {1, -Fe}; and b) the number j G {0, A^} of energy units 5 stored in the sensor ESD. 

States in which the sensor is idle, referred to as idle states, are instead denoted by Ij and they 
are uniquely characterized by the number j G {0, ...,A^} of energy units stored in the sensor 
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Figure 3. a) Discrete Markov chain used to model the evolution of the energy stored in the discrete BSD of a sensor in terms 
of the energy unit 5. In b. 1) and b.2) there are two outcomes of possible state transition chains for eg = 3. Grey shaded states 
indicate energy shortage condition. Some transitions are not depicted to simplify representation, (a — l~a andpc.fe = 1— Pc.fe)- 

ESD. EH is then associated to idle states given the assumption that any energy arrival in the 
current IR can only be used in the next IR (see Sec. III-C|) . 



B. Discrete Markov Chain (DMC) Model 

Operations of a sensor across IRs are as follows. When sensor Sm is not involved in an IR, it 
is in an idle state, say Ij, waiting for the next IR. When a new IR begins, the energy harvested 
in the last interval Tint is added, so that, if the ESD is not in energy shortage, the state makes a 
transition Ij — )• A} toward an active state, with I > es > j. Otherwise, if it is in energy shortage, 
it makes a transition Ij — )• // toward another idle state, with j < I < £5. If sensor Sm is not in 
energy shortage, it remains in state at the beginning of the IR only if it has a new measure to 
transmit, which happens with probability a. Instead, with probability a = 1 — a the state makes a 
transition toward an idle state as — )■ Ij. If there is a new measure, the sensor keeps transmitting 
it in successive frames until either the packet is correctly delivered to the EC, or its ESD falls 
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in energy shortage, or both. A collision in frame k happens with probability pc,fc = 1 — Pc,k (see 
Sec. IIV-A3I) and leads to a transition either Aj — )■ A^^l^, for j > 2es (no shortage after collision) 
or — )• /j-e^, for j < 285 (shortage after collision). Successful transmissions in frame k, which 
happens with probability Pcfc, instead leads to a transition — > Ij^es- Transition probabilities 
are summarized in Fig. IH where we have defined qj^N = Pr[-E//,m > iN—j)6] = 1 — Xlilo^"^ 
Note that, the probability a of having a new measure is only accounted for in active states in 
the first frame (i.e., in states Aj, for j G {0, A^}, see Fig. |4l-b)). In fact, being in any state Aj 
for k > 1 already implies that a new measure was available at the beginning of the IR. Notice 
that, according to the model above, state transitions in the DMC at hand are event-driven and 
do not happen at fixed time intervals. A sketch of the considered DMC is shown in Fig. [3]-a), 
while we show two outcomes of possible state transition chains in Fig. [31-b.l) and[3]-b.2). 

From Fig. |3]-a), it can be seen that, when go > 0, gi > and pc,k > 0, for k G {1, ...,Fs}, 
the DMC at hand is irreducible and aperiodic and thus, by definition, ergodic (see (TS\). In fact, 
if qi > 0, any state of the Markov model can be reached from any other state with non-zero 
probability, and therefore the Markov chain is irreducible. Moreover, the probability of having a 
self-transition from state Iq to itself is go > 0, and therefore state Jq is aperiodic. The presence of 
an aperiodic state in a finite state irreducible Markov chain is enough to conclude that the chain is 
aperiodic ifTSl Ch. 4, Th. 1]. Since the DMC is ergodic it admits a unique steady-state probability 
distribution = [0/^, ...,0/^^,0^!^, ...,0^Fe], regardless of the initial distribution, which can be 
calculated by resorting to conventional techniques (IS). This also guarantees the existence of 
limits © and ©. Vector represents the steady-state distribution in any discrete time instant 
of the interrogation period (i.e., during either any frames of an IR or idle period). However, 
to calculate Q and ^ we need the DMC steady- state distribution 0+ conditioned on being at 
the beginning of the IR. This can be calculated by recalling that between the end of the last 
issued IR and the beginning of a new one, sensor 5*^ can only be in any idle states Ij, with 
j G {0, ...,A^}, and thus its state conditional distribution 0~= [07,.. .,07 ,0^i ^---j^TfJj is 
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given hy (j)^, = (pi./ J2i=o 't'h^ ^ {0, N} and </) = 0, for all j, k. The desired distribution 
0+ of the state at the beginning of the next IR can be obtained as 0+ = 0^P, where P is 
the DMC probability transition matrix of the DMC in Fig. [3l-a) that can be obtained through 
Fig. m Note that, according to the transition probabilities in Fig. |4l starting from any state Ij, 
with j G {0, ...,N}, only states Ij, with j G {0, ...,es — 1} and states A], with j G {ss-, ■■■■,N} 
can be reached. Therefore, the only possibly non-zero entries of distribution 0+ are (jyj, for j G 



Once the DMC steady-state distribution 0+ at the beginning of any (steady-state) IR is 
obtained, we can calculate the steady-state distribution PE{n^oo){-) of the energy stored in the 
sensor ESD at the beginning of any (steady-state) IR, denoted by he = [7rE(0), ...,t:e{N)], by 
mapping the DMC states into the energy level set £ as follows 



The ccdf Ge( n->oo)(') is immediately derived from tce- Finally, we remark that analysis of FA 
and TDMA can be obtained by limiting the set of active states to A^^, A]y (i.e., no 

retransmission), and recalling that sensor Sm after transmission returns to idle states regardless 
of the success of transmission. 



{0, £5 - 1} and (l)\ifox j G {es, N}. 




(20) 
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Figure 4. State transition probabilities for the DMC model in Sec. IV-BI a) transition probabilities due to energy harvesting; b) 
transition probabilities due to the bidirectional communication with the FC. The transition matrix P can be derived according 
to the probabilities in a) and b) for all the values of G {1, F^} and j £ {0, A^}. 



VI. Backlog Estimation 

In this section we propose a backlog estimation algorithm for the DFA protocol (extension to 
FA is straightforward). Unlike previous work on the subject lfT6]||fT9l . here backlog estimation is 
designed by accounting for the interplay of EH, capture effect and multiple access. Computational 
complexity of optimal estimators is generally intractable for a large number of sensors even for 
conventional systems (see e.g., [fT9ll '). We thus propose a low-complexity two-steps backlog 
estimation algorithm that, neglecting the IR index, operates in every IR as follows: i) the FC 
estimates the initial backlog size Bi based on the ccdf Ge {s) of the ESD energy at the beginning 
of the current IR; ii) the backlog estimates for the next frames are updated based on the channel 
outcomes and the residual ESD energy. 

For the first frame, the backlog size estimate and the frame length are Bi = MuGe {b) and 



of Lk 



pBi 



, respectively. For subsequent frames, let us assume that the FC announced a frame 
pBk slots. The FC estimates the backlog size for frame A; + 1 by counting the number 
of slots that are successful (ND k) and collided {Nc,k) within the kth frame of length slots. 
Since the FC cannot discern exactly how many sensors transmitted in each successful slot, the 
estimate of the total number Cdm of sensors that collided in A^D,fc successful slots is CD,k = 
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{(^DM — 1) ^D,k, with /^Djt being the conditional average number of sensors that transmit in a slot 
given that the slot is successful (with no capture /3d, fc = !)• Similarly, for the collided slots we 
obtain Cc,k = (3c,kNc,k, where (3c,k is now conditioned on observing a collided slot. Derivations 
of (3r),k and /3c,k are in Appendix |Al Since the estimate of the total number of sensors that 
unsuccessfully transmitted is Ck = Cc,k + CD,k, the backlog size estimate Bk+i for the {k + l)th 
frame is obtained by accounting for the fraction of sensors within that are not in energy 
shortage: Bk+i = CkGE{{k + l)e\ke), where GE{{k + l)e\ke) = Pr [E„, > (A; + l)e\Em > ke]. 
The proposed backlog estimation scheme thus works as follows: 

{MaGE (e) if = 1 

(21) 
Ck^iGE{ke\{k~l)e) if k > 1 

Algorithm (|2T)) can be applied to any nth IR by deriving the ESD distribution PE{n)i-) (or 
GE{n)i-)) from any initial distribution pE{i)i-), by exploiting the DMC model in Sec. IV-BI 

VII. Numerical Results 

In this section, we present extensive numerical results to get insight into the MAC protocols 
design. Moreover, to validate the analysis proposed in Sec. |IV] and Sec. |Vl we compare the 
analytical results therein with a simulated system that does not rely on simplifying assumptions 
A.l and A.2. The performances of the backlog estimation algorithm proposed in Sec. IVll are 
also assessed through a comparison with the ideal case of perfectly known backlog at the FC. 

A. MAC Performance Metrics Trade-offs 

The energy EH,m{n) harvested between two successive IRs is assumed as geometrically- 
distributed with qi = Y'i[EH,m{n) = iS\ = ^(1 — ^)\ where ^ = + fin)- The average 
harvested energy normalized by e, referred to as harvesting rate, is E[EH,m{n)/ e] = ^h- 

The asymptotic time efficiencies dS]) for TDMA, FA and DFA protocols, are shown in Fig. [5] 
versus design parameter p (recall System performance is evaluated by considering: ijlh € 
{0.15, 0.35}, M = 400, 'jth = 'idB, a = 0.3; e is normalized to unity, energy unit is 5 = 1/50 so 
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that Es = 50 and = 10. We compare the analytical performance metrics derived in Sec. HVl with 
simulated scenarios for both known and estimated backlog. While the performance of TDMA is 
clearly independent of p, in FA and DFA there is a time efficiency-maximizing p, which is close 
to one (in |[8l the optimal value was p = I since the capture effect was not considered). The 
effect of decreasing (or increasing) the harvesting rate fin on the TDMA time efficiency is due 
to the larger (or smaller) number of sensors that are in energy shortage and whose slots are not 
used, while it is negligible for FA and DFA due to their ability to dynamically adjust the frame 
size according to backlog estimates Bk. The tight match between analytical and simulated results 
also validates assumptions A.l and A.2 and the efficacy of the backlog estimation algorithm. 

The asymptotic delivery probability, for harvesting rate fin G {0.05, 0.15, 0.35}, versus design 
parameter p is shown in Fig. |6] with the same system parameters as for Fig. [51 Unlike for the 
time efficiency, TDMA always outperforms FA and DFA in terms of delivery probability. In fact, 
sensors operating with TDMA and FA have the same energy consumption since they transmit 
at most once per IR, while possibly more than once in DFA. However, TDMA does not suffer 
collisions and thus it is able to eventually deliver more packets to the FC. The delivery probability 
strongly depends on the harvesting rate fin, which influences the ESD energy distribution and 
consequently the energy shortage probability. Moreover, DFA outperforms FA thanks to the 
retransmission capability when the harvesting rate is relatively high (e.g., fiH= 0.35), while for 
low harvesting rate (e.g., /i/fG {0.05,0.15}) DFA and FA perform similarly. In fact, for low 
harvesting rates, most of the sensors are either in energy shortage or have very low energy in 
their ESDs. Hence, most of the sensors that are not in energy shortage are likely to have only one 
chance to transmit, and thus retransmission opportunities provided by DFA are not leveraged. 

The trade-off between asymptotic delivery probability ^ and asymptotic time efficiency ^ 
is shown in Fig. |7] for different values of the harvesting rate fin G {0.05,0.15,0.35}. System 
parameters are the same as for Fig. [51 For TDMA, the trade-off consists of a single point on the 
plane, whereas FA and DFA allow for more flexibility via the selection of parameter p. When 
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increasing p more sensors might eventually report their measures to the FC, thus increasing the 
delivery probability to the cost of lowering time efficiency (see Fig. [5] and [61). For FA and DFA, 
the trade-off curves are obtained as maxp {pj[^^, s.t. pf^ = A for each achievable A. 

The impact of the capture effect on the performance metrics trade-offs is shown in Fig. [8l 
where we vary the SIR threshold ^th ^ {0.01, 3, 10}dB and keep the harvesting rate hh = 0.15 
fixed (other parameters are as in Fig. [5]). As expected, the lower the SIR threshold 'jth the higher 
the probability that the SIR of any of the colliding sensors is above 'jth, and thus the higher the 
performance obtained with ALOHA-based protocols. TDM A is insensitive to 'jth- 
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Figure 5. Asymptotic time efficiency ^ versus p, for different energy harvesting rates fiH G {0.15, 0.35}. Comparisons are 
between analytical derivations and simulated results with both known (Bt) and estimated backlog (Bt, see i2l\ ). (M = 400, 
7th = MB, a = 0.3, = 10, e = 1, 5 = 1/50). 
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Figure 6. Asymptotic delivery probability (O versus p, for different energy harvesting rate pn £ {0.05,0.15,0.35}. 
Comparisons are between analytical derivations and simulated results with both known (Bk) and estimated backlog (Bk, see 
inj), (M = 400, 7th = SdB, a = 0.3, F, = 10, e = 1, 5 = 1/50). 
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VIII. Conclusions 

The design of medium access control (MAC) protocols for single-hop wireless sensor networks 
(WSNs) with energy-harvesting (EH) devices offers new challenges as compared to the standard 
scenario with battery-powered (BP) nodes. New performance criteria are called for, along with 
new design solutions. This paper addresses these issues by investigating the novel trade-off 
between the delivery probability, which measures the capability of a MAC protocol to deliver 
the measure of any sensor in the network to the intended destination (i.e., fusion center, FC) and 
the time efficiency, which measures the data collection rate at the FC. The analysis is focused on 
standard MAC protocols, such as TDM A, Framed- ALOHA (FA) and Dynamic-FA (DFA). Novel 
design issues are also discussed, such as backlog estimation and frame length selection. Extensive 
numerical results and discussions validate the proposed analytical framework and provide insight 
into the design of EH- WSNs. 

Appendix A 

Average Number of Sensor Transmissions per Time-slot 

The conditional averages (3£),k and (3c,k are calculated similarly to |[8l by accounting for the 
capture effect and an arbitrary p. Let Y be the number of simultaneous transmissions in the 
same slot, and let Uk and Ck respectively be the event of successful and collided slot in frame 
k, the average number of sensors per successful and collided slot are respectively 

oo oo 

f3D,k = J] J Pr [F = J \Uk] ; ^c,k = J] j Pr [F = j\Ck] (22) 

To calculate ^ consider A.l and A.2 and allow the number of possible interfering users up to 
infinity as in Sec. lIV-A2l By exploiting the Bayes rule, we have Pr [Y = j\Uk] = Pr[Uk\Y = j] 
where Pr [l(k\Y = j] = jpc,k{j - 1), Pr [Y = j] = e~p/{p^j\) and Pr [Uk] = p?J^ (see[H. We 
can similarly obtain I3c,k given that Pr [Ck] = 1 — Pr [W^] — /3 (0, B, L), where /3 (0, B,L) e p 
is the probability of an empty slot, and Pr [Ck\Y = j] = \ — Pr \Uk\Y = j] for j > I. 
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